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ABSTRACT 



The direct current resonance probe technique was applied to a 
moderately high density positive column xenon discharge in the presence 
of moving striations. The resonance frequencies obtained were of the 
order of 500 megacycles per second, corresponding to electron plasma 
frequencies of the order of one to one and a half gigacycles per second. 
The results of the experiment indicate that the resonance probe tech- 
nique may be a useful method for the determination of electron tempera- 
ture and density, but not electron-neutral collision frequency, in 
moderately high density plasmas. Experimental difficulties are dis- 
cussed, and recommendations for the course of future work with 



resonance probes are presented. 
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1 . 



Introduction. 



The resonance probe technique was first developed in I960 by 
Takayama, Ikegami, and Miyazaki [1], as a means of quickly deter- 
mining electron temperature, density, and collision frequency in a 
plasma. 

The method consists basically of applying a high-frequency signal 
to a Langmuir probe and noting the behavior of the collected direct 
current. Referring to Figure 1, the direct current shows three dis- 
tinct regions as the frequency of the applied signal is varied: the first 

is an increase in the dc level in the region below the resonance fre- 
quency, independent of the signal frequency, from which the electron 
temperature can be determined; a noticeable peak at the resonance 
frequency; and at higher frequencies, a rapid fall-off of the direct 
current to the level corresponding to no applied signal. From the 
position of the resonance frequency, electron plasma frequency can be 
calculated, and from the half-width of the peak, the electron-neutral 
collision frequency can be determined directly. The basic circuit for 
the dc resonance probe method is illustrated in Figure 2. 

In this experiment, the resonance probe technique was applied to 

a moderately high density (on the order of 10^ electrons cm xenon 

-3 

positive column with pressures on the order of 10 torr) . 

There were major differences in the type of experimental conditions 
here as compared to most previous work. Most previous work has been 
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done in "back diffusion" plasmas, with electron densities on the order 
6 -3 

of 10 cm , and the type of gas has usually been mercury or cesium 
vapor. In the present investigations, in addition to the high density 
xenon, there were moving striations present in all ranges covered. These 
striations had wavelengths of seven or eight cm and frequencies of 
25 - 100 kilocycles per second, with no apparent linear dependence on 
plasma parameters. They are ionization type waves which usually move 
from the anode to the cathode through the positive column [13]. 
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2. Resonance Probe Theory. 

Takayama and Ikegami [2] believed that the resonance frequency 
occurred at the electron plasma frequency, which is a function of electron 
density only: 

(1) U) = ( 4IIn l e l \ = 2n • 8.97xl0 3 (n(cm -3 )) 1/2 , 

p V m y 

where uu is the electron plasma frequency, n the electron density (not 
P 

necessarily constant), and e and m the electron charge and mass, 
respectively. Subsequent investigations [5] [7] [8] [9], however, have 
shown that the dc resonance peak always occurs below a) 

For superposed signals of form A* sin ujt , below the resonance 
frequency, the increase in direct current is explained as follows. The 
electron current density, , to the Langmuir probe, is given by the 
well-known relation: 



e V 




where K is the Boltzmann constant, T the electron temperature, and 
V is the potential difference between the plasma and the probe. The 
subscript ^ will refer to the current density due to the probe bias, V, 
which is held constant in resonance probe investigations. 

The electron velocity distribution function is also assumed to be 
isotropic, and is determined by: 
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(3) 



f(v) 



V 



4m 

T 

e 




where v is the electron velocity. 

According to Ikegami and Takayama [2], the direct current increase, 
6j , is then given by the second derivative of the electron current density 
versus potential difference curve: 



(4) 



1 A 2 

6J = 4 A 2 
dV 



If the distribution is assumed Maxwellian, the dc increase can be 
written: 




where I (z) is the zero order modified Bessel function, and A is the 
o 

amplitude of the superposed signal. From this the electron temperature 
can be determined: 




For frequency ranges above the resonance frequency, the decrease 

in 6j to zero is explained by Ichikawa [3] as being due to the effect of 

the oscillating potential cancelling out as an average, since the electrons 

cannot follow the fluctuating electric field, the transit time, t , of the 

o 

electrons across the probe sheath, is much longer than the period of 
the potential oscillation, o> ^ . 

Solving the equation of motion for the electrons: 
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( 7 ) 



mx = e E + e E sin (ud t - a ) , 
o 1 



V A 

where E = — , E = — , s is the sheath thickness and a is a phase 
o s 1 s 

angle due to the plasma impedance defined by tan a = e / e , , where 
€ and €. are the real and imaginary components of the plasma permit- 
tivity, respectively; and x is the distance of the electrons from the probe. 
Considering an electron initially at the edge of the sheath, which will 
require a minimum velocity v^ (hereafter called the "critical velocity") 
to reach the probe, leads to the following initial and boundary conditions 
for the solution of the equation of motion: 

x = s, x = - v at t = 0, and x = 0 at x = 0 . 



Then the critical velocity is derived: 



e E 

(8) v ( a , u> ) = - t 

c mo 



e E, 



(o) + fcos a - cos(o) t (a) - a ) ") . 

m a) V o y 



Then: 



(9) 



} 2n 

j o + 6j = 2 J J v^ (a , a) ) vf(v) dvda 



Solving equation (8) for the critical velocity in the low frequency 
region yields: 



( 10 ) 



<»> ■ ci < 



V - A sin cc 



o 



1/2 



Substitution of equation (10) into equation (8) leads to equation (5). 
Integrating equation (7), one obtains: 
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s - 



( 11 ) 



2 



e E 

o 

m 



2 

t 



e E 



mo) 



cos (a) t - a ) 



e E 



mu) 



(^sin a 



+ 



sin (cut - oi 




0 . 



In the high frequency limit, the transit time of electrons across the 
sheath is found from equation (11): 



( 12 ) 



T 

O 



lim (t) = ^ 

it) — ^ » 



2ms 
e E 

o 



) 



1/2 



Substituting equation (12) into equation (8), yields: 




from which 6j approaches zero as a) s becomes greater than v. 

In the frequency region at and near the resonance frequency, the 
mechanism is not as straightforward, and no universally accepted theory 
is available at present. A summary of some of the existing theories 
which are applicable follows. 

The plasma is assumed to be homogeneous, and for simplicity in 
development, plane parallel probe geometry is used. The plasma 
system can be represented schematically as in Figure 3 (a) , with the 
high frequency potential applied across the probe and a reference elect- 
rode, usually the cathode. The resulting equivalent circuit is illus- 
trated in Figure 3(b), after Levitskii and Shashurin [4], Uramoto, et al 
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[5], and Mayer [6]. 

It has been confirmed that both the rf and dc current components 
show a maximum below the electron plasma frequency, and the the rf 
current drops to minimum near the electron plasma frequency. The 
former case (maximum) is interpreted as a series resonance, and the 
latter as a parallel resonance. Uramoto, et al [5], showed that the 
rf component increased monotonically at frequencies above the electron 
plasma frequency, while the dc component, 6j, as previously discussed, 
drops rapidly to zero. 

As the present investigations were concerned primarily with the 
dc resonance, the development will follow Ichikawa [3] and Ikegami 
and Takayama [2], realizing that the theory must be applied to the 
probe-plasma series capacitive circuit, and that the resonance fre- 
quency, id , is not the same as the electron plasma frequency, 

If an external oscillating potential is applied at a frequency near 

the electron plasma frequency, the much slower ions present in the 

sheath about the biased probe will not be able to follow the oscillating 

electric field and their distribution will be a function only of the probe ^ 

bias, V. The electrons in the sheath, however, will follow the applied 

oscillating electric field and the resulting charge fluctuation in the 

sheath acts as an oscillating dielectric layer. The external perturbing 

field, E , can be expressed as: 
ext 

A 

(14) E ext (x, t) = — D(x) sin (uot - a ) , 
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where 



(15) D(x) - i [ 



CO sinks ikx 

— ; e dk , 



where k is the wave number. When these latter equations are used to 
solve the linearized Boltzmann equation: 



(16) || + v|i ♦ i 

at ax m a v 



— E + — - + 

m ext 3 v at 



coll , 



where f is the distribution function in an equilibrium state and f is 
o 

the perturbation due to the applied signal, and also 



(17) 



aE 

ax 



= 4 



n e J 



f dv , 



the following solution is obtained: 



af 



(18) f (x, v, t) = — — - A f duo e * ^ ° t f o' ( u) +oo)e 

mavj o \ o 



- a 



-a(oJ o -cjo)e v 



1 

2n J 



dk e 



ikx 



sin k s 



oo-kv+iv € (k, a) ) ks 

o c o 



where a> is taken as a dummy variable, and v is the effective electron 
o c 

collision frequency. 

When the mean path is much greater than the sheath thickness, the 
integral over k in the last equation can be estimated asymptotically: 



(19) Lim 



if i 



dke 



ikx 



1 



1 sin k s 



1 



1 1 



o>-kv+iv e (k , a) ) ks o> + iv e(o,u) ) s 

CO CO 



s o 
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Due to the time dependence of the perturbing rf signal, equation (9) is 
evaluated at t = 0, and assuming Maxwellian distribution of velocities, 
the total distribution function becomes: 



{20) f total (V) = f o (v) + f(x ’ v; t = 0) 



1 lA 

2 KT 



v 



f 

o 



(v) 



1 

OJ s 



2 

(D 



2 

(U + 

p 



V 



2 

c 







cos a 




The real and imaginary components of the dielectric permittivity are 
defined by: „ 



(U 

(21a) ® r (u>) = 1 - 2 

0 ) - v 

c 



and 



(21b) 







2 tu v 

c 



2 2 ' 
U) - V 

c 



Substituting equation (20) into equation (9) yields the final result: 

eV eV 

,,,, . KT , , , , . KT “p^D u> 

(22) j o + 6j(u.) = J c e I 0 (») + J 0 « — j 2 

(JU + V 



B , 



where: 
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B 





v e 
c r 



CD 

O 




+ ) 




+ 







where X ^ is the Debye length, and: 



(23) 



T / \ 1 r2n - 1/2 X c 

X o <UJ) = 2Ti J e da 



1/2 

where the dimensionless critical velocity, x (cv, u)) = (m/KT) v (a, cd), 

c c 

and 



/" S M N 


\ - _L_ 


r 2n 


r oo 


V c(u)) ^ 


f 2n 


J o ■ 


J 

X 



2 - l/2x 

x e 



^sin Qf\ 
dx vda 

V COS Ot J 



Equation (22) is a general solution for total dc collected by the 
probe. Using the critical velocity for the low frequency limit given by 
equation (10), the first term of equation (22) reduces to equation (5). 

If the approximation 



(25) 




x e 
c 



1/2 x c 2 



is taken, equation (22) reduces to: 



(26) 6 j 



2 




2 v / a) 
c 

2 



0) 






e A 
KT 
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Equation (26) can also be derived by solving the fundamental equations 
for linear response of the plasma electrons, equations (16) and (17) , 



directly for dc current density, assuming that 



dt coll 



- V f , 

c 



eV/KT is much greater than one, that A/V is much less than one, and 
2 

that ( v / cu ) is also much less than one. 
c 

Furthermore, if equation (26) is solved for the resonant peak 

height, 6j , assuming a) = cd = cu : 
r r p 




where I (z) is the first order modified Bessel function. 

The half width of the resonance peak determines the effective 
electron-neutral collision frequency, the relation as shown by 
Takayama, et al [2], is: 

(28) A (D = 2 v 

c 

This has been experimentally verified to within an order of magnitude 
by Cairns [7], using a low density mercury vapor plasma. 
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3. Relationship of Resonance Frequency to Electron Plasma Frequency. 



As previously noted, the resonance peak of rectified current always 
occurs below the electron plasma frequency, and that this is due to 
series representation of the sheath-plasma interaction. 

Mayer [6] developed a theory for plane probe geometry and cal- 



culates the complex permittivity 6 , of the plasma capacitor as: 

e xx 

2 2 

2 2 u)-o) + i v u> 

< 29 > * eff * ‘r + *1 * 1 + <% - «„ > — V 2 2 



(( 1 ) - 0 ) ) 4- v w 

r c 



where 

(30) 



2 2s 2 

u) = 

r p + 2s p 



This model has been experimentally confirmed by von Gierke, et al [8] 
for a low density cesium vapor plasma with unequal sheath thicknesses 
at the probe and reference electrode ( s and s , respectively) . 

X £* 

Harp and Crawford [9] expanded this theory to spherical geometry 
and calculated: 



(31) 




( 1 + 



R 



a X 



D 



2 

0) 

P 



where R is the spherical probe radius and a X replaces s, and a is 
a constant of proportionality. 
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4. Experimental Results and Conclusions. 

The resonance probe technique was applied to a xenon positive 
column discharge in a 60 cm Pyrex tube of 0. 8 cm inside diameter. 

The probe was a 0. 195 mm radius right-angle cylindrical tungsten probe 
located halfway between an oxide coated cathode and a tungsten rod 
anode. Neutral gas pressures of one to four microns of mercury were 
used in most of the experiments with electron densities in the plasma of 
the order of 10^ cm ^ . Tube currents were 20, 30, and 50 milli- 
amper es. 

Early attempts to find a consistent resonant increase in collected 
direct current were unsuccessful. It was believed that, in accordance 
with Harp and Crawford [9], the probe radius was too small in relation 
to \ ^ , thereby damping out any resonant peak. In order to eliminate 
this problem, a second probe was installed in the system near the right - 
angle probe, towards the anode. This was a thin tantalum band, two mm 
wide, which was mounted around the inner wall of the tube. The 
apparatus is diagrammed in Figure 4. After installation of this band, 
resonance peaks were observed, both on it, and on the original cylindrical 
probe, when either probe was at floating potential and the other biased 
near floating potential and with an applied rf signal. It would appear 
that the presence of the second conducting surface is essential to ob- 
servance of resonance frequencies in this type of discharge in that it 
provides a closer reference electrode than either the anode or cathode. 
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The radial profile of the electron density is assumed to be parabolic, 
and that the electron density at the wall is approximately 60% of that at 
the center of the tube. A typical plot of 6 j/j Q versus the applied 
frequency is shown in Figure 5, with results obtained from it listed in 
Table 1. 

It was expected that density variations caused by moving striations 
would broaden the resonance peak. This was confirmed in all experi- 
mental conditions as illustrated in Table 2. For calculation of v , the 

c 

electron temperature as obtained from d(ln j^) / dV plots was inserted 
into 



(32) 



u /8KT \l/2 

= — = (n^J (p c P> 



where p is the neutral gas pressure, in torr, and is the probability 
of collision at one torr, [10], u is the average electron velocity, and 
X ^ is the electron mean free path. 

The expected peak height, 6 j / j , as calculated from equation (27), 
did not show a consistent relation to those measured, but, instead, 
varied randomly. It was confirmed, however, that the magnitude of the 
peak height, and the dc increase observed below the resonance frequency, 
were proportional to the amplitude of the applied signal. The results 
of peak height comparisons are given in Table 3, the amplitude of the 
applied signal was 0. 7 volts. 

The results of the measurement of the resonance frequency with 
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respect to the calculated electron plasma frequency appear to be closer 

to the model of Mayer [6] than to that of Crawford and Harp [9]. As 

shown by von Gierke, et al [8], the sheath thickness is assumed to be 

a function only of the Debye length, A. , i. e. > s = a\^ . 

For the right-angle cylindrical probe (CP), the value of a is 2. 66 

for closest fit of the ratio of measured uo to calculated uo , using 

r r 

2 

equation (30) and u> obtained from d(j Q )/dV plots. For the wall, or 

ring, probe (RP), the value of a is 3. 5. In the latter case, uo is 

P 

approximated by determination of the random ion current density, j . , 
from Langmuir probe curves, and the relation from Glas stone and 
Lovber g [11]: 



(33) 



n. 

l 



r l 



( 0. 40 




1/2 



where m is the ion mass. This assumes charge neutrality, i. e. , the 
ion density, n. = n. 

A sample of results of the above measurements is given in Table 4, 

and the ratio of calculated uo to measured uo is plotted against experi- 

r r 

mentally measured uo in Figure 6. 

In the frequency region below the resonance frequency very poor 
agreement with theory was obtained. The increase in the dc level ap- 
peared to be independent of applied signal frequency, but a large random 
scatter was noticed, and, at best, only an estimate of 6 j / j could be 
made. This was most probably due to drift of j and to the difficulty in 
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accurately reading the changes in the total current density, j , since 
6 j = j ~ J Q was very small. Sample values are listed in Table 5. 

In the region above the resonance frequency, 6 j did not drop rapidly 
to zero, and appeared to have one or more resonance peaks. These 
peaks, as suggested by Ikezi and Takayama [12], may be Tonks - 
Dattner resonances. 

Since a resonant peak can be detected even at these comparatively 
high densities, it is concluded that the method of resonance probes may 
be useful as a diagnostic technique in positive column discharges for 
determination of electron temperatures and densities. However, due to 
the presence of moving striations, it is believed that the method may 
yield an average collision frequency of the plasma. 
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5. Recommendations for Future Investigations. 

In future work it is recommended that the rf current be examined 
with the ultimate goal of developing a simple method of direct determi- 
nation of electron plasma frequency. 

It is further recommended that the dc resonance probe again be 
attempted in the absence of a second close electrode. The problem of 
electron temperature determination using the resonance probe technique 
may be overcome if a leveled-power rf generator is obtained, and a 
method can be devised to eliminate the drift of j , and to measure only 
6 j rather than j . The apparent resonant peaks above oj should be more 
carefully investigated to determine whether or not they do actually 
represent Tonks - Dattner resonances or some other phenomenon. 



25 



BIBLIOGRAPHY 



1. Takayama, K. , H. Ike garni, and S. Miyazaki. Phys. Rev. Letters, 

5, p. 238, I960. 

2. Takayama, K. , and H. Ikegami. Resonance Probes , Institute of 

Plasma Physics, IPPJ-10, Nagoya University, Nagoya, 

Japan, March, 19 63. 

3. Ichikawa, Y. H. Theoretical Analysis of the Characteristics of 

Resonance Probe, Institute of Plasma Physics, IPPJ-14, 

Nagoya University, Nagoya, Japan, October, 1963. 

4. Levitskii, S. M. , and Shashurin. Sov. Phys. - Tech. Phys. 8, 

p. 319, 1963. 

5. Uramoto, J. , J. Fujita, H. Ikegami, and K. Takayama. RF Current 

Component in the Resonance Probe , Institute of Plasma Physics, 
IPPJ-19, Nagoya University, Nagoya, Japan, October, 1963. 

6. Mayer, H. M. Measurements with a Wide-Band Probe, Proceedings 

of the Vlth International Conference on Ionization Phenomena in 
Gases (Paris, July 1963) , Vol. 4, 1964, pp. 129-134. 

7. Cairns, R. B. Measurements of Resonance Rectification using a 

Plasma Probe, Proc. Phys. Soc. , Vol. 82, 1963, pp. 243-251. 

8. vonGierke, G. , G. Muller, G. Peter, and H. H. Rabben. On the 

Influence of Ion Sheaths upon the Resonance Behaviour of a rf 
Plasma Probe, Zeitschrift fur Naturfor schung, Vol. 19, May, 
1964, 1107-1111 - 

9. Harp, R. S. , and F. W. Crawford. Characteristics of the Plasma 

Resonance Probe, Jounral of Applied Physics, Vol. 35, 
December, 1964, pp. 3436-3446. 

10. Brode, R. B. Revs. Modern Physics, 5, 1933, p. 257. 

11. Glasstone, S. and R. H. Lovberg. Controlled Thermonuclear 

Reactions, D. van Norstrand Company, Inc., I960. 



26 



12 . 



Ikezi, H. , and K. Takayama. Resonances of Radio Frequency Probe 
in a Plasma, Institute of Plasma Physics, IPPJ-43, Nagoya 
University, Nagoya, Japan, March, 1966. 

13. Cooper, A. W. Experiments on the Origin of Moving Striations, 
Journal of Applied Physics, Vol. 35, 1964, p. 2877. 



27 




FIGURE 1 



Theoretical Plot of the Increase of dc Current 
Versus Frequency of the Applied rf Signal 
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FIGURE 2 



Basic Circuit for the Resonance Probe 
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FIGURE 3 



(a) The Model Proposed for the Plane Probe 
Plasma- Sheath System 

(b) The Corresponding Equivalent Circuit 
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(1) Signal Generators: (2) Ammeters: 

10-480 me Hewlett-Packard 608 c Weston Model 430 

450-1200mc Hewlett-Packard 612a DC Mi cro ammeter 

10kc-50 me SG-85 / URM - 25P Keithly Instruments 

Model 600 A Electrometer 
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Typical Plot of 6j/j Versus Signal Frequency Obtained by the Resonance Probe Technique. 

The gas is xenon, tube current 20 milliamperes, pressure 2. 15 microns 
of mercury, and the probe used was the axial cylindrical probe. 



TABLE 1 



Analysis of Figure 5. Basic parameters are: xenon positive 
column discharge, pressure of 2 . 15 microns of mercury, and tube 



current 20 milliamper es. 




Resonance frequency, a) ^ / 2 11 : 


490 megacycles per second 


Half width, A u) / 2 n : 


100 megacycles per second 


Collision frequency, v ^ / 2 II : 


50 megacycles per second 


Resonance peak height, 6 j ^ : 


0. 090 


Electron temperature, KT/e : 


5. 0 5 ± 1. 50 ev 



From the corresponding j versus V probe curve, the following compara- 



tive data is obtained: 
Electron density, n: 


2. 22 x 10^ cm 



Electron plasma frequency, uu / 2 : 1320 megacycles per second 

P 



Collision frequency: 


32, 1 megacycles per second 


Electron temperature : 


4. 05 ev 


Electron Debye length, 


0. 10 mm 


* 

Probe potential, V : 


23 volts 



33 




o 

o 

r- 



o 

in 

xQ 



O 

O 

v£> 



O 

lO 

in 



o 

o 

to 



o 

LO 



o 

o 



CM 




O 

i — i 

Ui 



cn 

CM 



0] 

CM 



a 

w 



a 



34 



Plot of a) ^ (Calculated) / u) (Measured) versus u) / 211 (Measured by 
Resonance Probe Technique, in Megacycles per Second), 
x - from ring probe data ; o - from cylindrical probe data . 



Representative Values of Measured Electron-Neutral 
Collision Frequencies and Comparison with Calculated Values 
(All frequencies are in megacycles per second. ) 
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TABLE 3 



Representative Values of Measured Resonance 
Peak Height and Comparison with Calculated Values 



Electron Density 




Peak Height ( 6 j / j ) 
r o 

Measured (1) Calculated (2) 



Ratio of 
Column 2 to 1 



Cylindrical Probe 



2. 22x 10 10 


0. 046 


0. 122 


2. 65 


2. 22x 10 10 


0. 090 


0. 105 


1. 17 


2. 18x 10 10 


0. 119 


0. 098 


0. 82 


3. 70xl0 10 


0. 078 


0. 158 


2. 03 


4. 36x 10 10 


0. 118 


0. 046 


0. 39 




Ring Probe 






1. 03xl0 10 


0. 038 


0. 0305 


0. 80 


1. 32x 10 10 


0. 081 


0. 0373 


0. 46 


2. 01 x 10 10 


0. 053 


0. 0732 


1. 38 


2. 57x 10 10 


0. 081 


0. 0425 


0. 525 
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TABLE 5 



Representative Values of Electron Temperature 
as Determined by the Low Frequency Signal DC Increase 
Compared with Corresponding Values from Langmuir Probe Curves 



Pressure 

Current (microns of 

(ma) mercury) 



Electron Temperatures 

Resonance Langmuir 

Probe Probe 



Cylindrical Probe 



20 


1. 6 


9. 82 


4. 30 


20 


2. 15 


5. 05 


4. 05 


20 


2. 4 


3. 53 


3. 80 


20 


4. 5 


3. 37 


3. 62 






Ring Probe 




20 


1. 4 


7. 85 


4. 74 


20 


3. 0 


3. 54 


4. 34 


20 


4. 5 


6. 43 


4. 34 


50 


2. 3 


3.93 


5. 63 
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